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P  REFACE 


This  Memorandum,  the  first  in  a  series,  presents  an  adaptive 
modal  for  scheduling  cargo  flights  within  the  Military  Airlift  Command 
(MAC).  It  provides  the  technical  details  for  a  pair  of  mathematical 
models  developed  for  this  scheduling  problem.  The  models  employ  an 
adaptive  technique  that  copes  with  the  uncertainty  associated  with 
some  of  MAC's  requirements. 

The  study  arose  from  a  joint  RAND/MAC  interest  in  dealing  with 
the  volatility  that  appears  in  MAC's  flight  structure,  particularly 
that  resulting  from  special  assignment  airlift  missions. 

This  Memorandum  is  primarily  addressed  to  operations  analysts  and 
schedulers  in  MAC,  but  should  also  interest  those  concerned  with  long- 
run  planning  of  airlift  routing  and  those  responsible  for  procuring 
commercial  lift  both  at  MAC  and  in  the  DOD.  Future  Memorandums  will 
provide  a  nontechnical,  user-oriented  approach  and  document  the  computer 
programs  for  the  algorithm  described  here. 

A  description  of  the  approach  developed  here  was  briefed  to  the 
Commander  of  MAC  and  his  staff  in  February  1968. 
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SUMMARY 


This  Memorandum  describes  an  analytic  approach  to  flight  scheduling 
within  the  Military  Airlift  Command.  The  model  takes  explicit  account 
of  the  uncertainty  present  in  the  cargo  requirements  or  demands  con¬ 
fronting  MAC.  For  computational  feasibility,  the  approach  consists  of 
two  related  models:  (1)  a  monthly  planning  model  that  produces  an 
initial  schedule  and  (2)  a  daily  model  for  making  periodic  changes  in 
the  schedule.  They  are  intended  to  schedule  both  the  21st  and  22d 
Air  Forces  together  and  take  account  of  interactions  between  them. 

Both  are  formulated  as  two-stage  stochastic  linear  programs.  A  de¬ 
tailed  mathematical  description  of  each  model  and  its  physical  inter¬ 
pretation  is  given. 

The  monthly  model  determines  the  number  of  flights  —  channel  and 
special  assignment  airlift  for  each  type  of  aircraft  in  the  MAC  fleet. 
Excess  demands  on  certain  routes  are  assumed  to  be  met,  at  least  in 
part,  by  additional  spot  procurement  of  commercial  lift.  The  flight 
assignment  is  determined  by  minimizing  the  expected  total  system  cost, 
which  consists  of  operating  costs,  costs  of  reallocating  aircraft  to 
different  routes,  commercial  procurement  costs,  and  other  penalty  costs 
of  excess  demand.  The  model  accounts  for  limitations  on  the  number  of 
flying  hours  and  the  carrying  capacities  of  various  aircraft  in  satis¬ 
fying  demands. 

In  the  daily  model  the  number  of  aircraft  of  each  type  to  svitcli 
from  one  route  to  another  and  the  number  of  commercial  flights  on  spot 
contract  to  add  on  the  current  day  are  the  principal  decision  variables. 
These  are  determined  by  balancing  operating,  procurement  and  redistri¬ 
bution  costs  against  the  expected  costs  of  additional  cargo  delay.  The 
current  state  of  the  cystem  --  the  amount  of  unmoved  cargo  on  various 
routes  and  the  position  of  aircraft  throughout  the  system  --  plays  a 
role  in  determining  these  decisions. 

These  scheduling  models  represent  the  first  known  attempt  to  im¬ 
plement  the  general  methodology  of  stochastic  linear  programming  on  a 
large-scale  operating  problem.  Consequently,  a  description  of  two  vari¬ 
ants  of  an  algorithm  recently  developed  for  this  class  of  problems  is 
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presented.  Both  versions,  which  use  ideas  from  convex  programming, 
make  extensive  use  of  linear  programming  codes  for  the  brunt  of  the 
calculations.  The  models  may  thus  be  solved  by  augmenting  existing 
linear  programming  routines.  A  computer  program  for  solving  a 
large  system  is  currently  under  development. 

Although  the  model  is  specifically  designed  for  the  MAC  system, 
it  is  hoped  that  with  appropriate  modifications  it  may  be  applicable 
to  other  scheduling  problems  in  which  uncertainty  is  a  factor. 
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I .  INTRODUCTION 

This  Memorandum  provides  the  technical  details  of  a  mathematical 
model  for  scheduling  flights  in  the  Military  Airlift  Command  (MAC) . 

The  model  takes  explicit  account  of  the  uncertainty  surrounding  the 

requirements  at  the  time  scheduling  decisions  must  be  made,  as  well  as  j 

of  the  various  physical  constraints  on  the  system.  The  proposed  models  j 

! 

are  adaptive  in  the  sense  that  current  information  about  the  state  of  j 

the  system  (cargo  levels  at  various  origins,  distribution  of  vehicles  • 

throughout  the  system)  is  used  periodically  in  making  decisions  regarding  1 


changes  in  the  flight  schedule.  1 

The  approach  discussed  here  consists  of  two  components  —  a  monthly 
and  a  daily  flight  planning  model.  Both  are  formulated  as  stochastic 
two-stage  linear  programs,  and  can  be  solved  by  algorithms  based  upon  j 

standard  linear  programming  calculations.  Consequently,  the  model  can  | 

be  applied  to  large-scale  systems.  This  should  permit  scheduling  the  \ 

i 

MAC  system  as  a  whole  rather  than  the  21st  and  22d  Air  Forces  separately,  * 

with  transfers  of  aircraft  between  them  on  an  ad  hoc  basis,  as  is  cur-  J 


rently  the  case. 

The  Memorandum  is  organized  as  follows.  The  next  section  contains 
a  brief  description  of  the  scheduling  problem  facing  MAC  and  discusses 
the  overall  strategy  employed  in  treating  it  analytically.  This  is 
followed  in  Sec.  Ill  by  a  mathematical  statement  of  the  monthl*  and 
daily  models,  which  are  similar  in  genei *1  structure.  Section  IV  con¬ 
tains  a  discussion  of  two  variants  of  an  algorithm  for  solving  these 
models,  preceded  by  a  short  discussion  of  pertinent  concepts  for  two- 
stage  linear  programming  under  uncertainty.  In  Sec.  V  we  present  a 
small  numerical  example  for  the  monthly  model. 
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II.  PROBLEM  STATEMTjjT 

DESCRIPTION  OF  MAC  SYSTEM 

Airlitc  missions  within  MAC  are  of  two  basic  types:  channe I 
missions  and  special  assignment  airlift  missions  (SAAMS) .  Channel 
missions  fly  on  a  regular  basis  between  fixed  origins  and  destinations: 
MAC  bases  in  the  continental  United  States,  aerial  ports  cf  embarkation 
(APOEs)  and  overseas  bases,  and  aerial  ports  cf  debarkation  (APODs) . 

For  example,  flights  from  Travis  Air  Force  Base  to  Bangkok  are  channel 
missions.  They  carry  cargo  (on  a  priority  basis)  and  armed  services 
passengers.  They  are  scheduled  in  advance  for  each  calendar  month, 
although  the  schedules  (number  of  flights  and  departure  times)  are 
changed  during  the  month,  often  daily,  based  upon  deviations  from 
expected  cargo  generation.  Many  of  these  flights  make  intermediate 
stops  for  refueling  and  crew  changes.  Channel  cargo  is  shipped  into 
APOEs  for  storage  until  shipped. 

SAAMS  are  flights  from  points  of  original  cargo  generation,  such 
as  factories  or  warehouses,  to  points  at  or  near  the  ultimate  consignee. 
Such  flights  arise  on  an  irregular  basis,  and  are  scheduled  if  sufficient 
advance  notice  is  provided.  The  time  and  place  of  cargo  pickup  is 
usually  specified.  Often  these  flights  arise  after  the  beginning  of 
the  month  and  require  a  temporary  diverting  of  aircraft  from  channel 
missions,  since  certain  SAAMS  are  designated  higher  priority  than  chan¬ 
nel  traffic.  The  users  of  a  SAAM,  e.g.,  a  given  military  organization, 
must  in  effect  charter  an  entire  aircraft,  although  MAC  decides  what 
type  of  aircraft  will  be  supplied. 

MAC  is  divided  into  two  air  forces,  the  21st  and  22d .  The  22d 
Air  Force  comprises  APOEs  on  the  West  Coast  and  serves  points  in  Asia. 

The  21st  comprises  APOEs  on  the  East  Coast  and  serves  bases  in  Europe, 
the  Middle  East,  and  South  America,  although  currently  several  channels 
serve  bases  in  Southeast  Asia.  MAC  Headquarters  assigns  SAAMS  to  the 
two  Air  Forces;  those  originating  west  (east)  of  the  Mississippi  River 
are  usually  assigned  to  the  22d  (2!st)  Air  Force.  Often  aircraft 
belonging  to  one  Air  Force  are  temporarily  assigned  to  the  other  Air 
Force  to  meet  these  demands  when  the  type  of  aircraft  required  is  in 
short  supply. 
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SCHEDUL1NG  PROBLEM 

In  carrying  out  its  operations,  MAC  is  faced  with  random  fluctu¬ 
ations  in  the  amount  of  cargo  and  passengers  to  be  carried  over  its 
various  routes. 

There  is  monthly  uncertainty  because  shippers'  forecasts  of  require¬ 
ments  for  the  month  as  reported  to  MAC  at  various  lead  times  differ 
significantly  from  the  actual  amount  of  requirements  shipped.  Secondly, 
the  amount  of  cargo  actually  received  during  a  day  is  partly  unpredict¬ 
able  due  to  random  variability  in  production  schedules  for  manufactured 
items,  delays  in  transit  and  depot  processing,  and  other  factors. 

A  scheduling  problem  arises,  since  demand  is  not  known  with  cer¬ 
tainty  and  aircraft  differ  with  respect  to  their  payloads,  flying  times 
and  operating  costs  over  different  routes.  Therefore,  aircraft  of  dif¬ 
ferent  types  must  be  assigned  to  various  routes  over  the  monthly  planning 
period.  In  addition,  MAC  can  obtain  additional  lift  from  commercial 
sources  on  call  contracts  to  handle  excessive  demands  on  some  routes. 
Furthermore,  the  assignments  of  aircraft  cannot  violate  constraints 
for  flying  hours. 

As  time  unfolds  and  demand  becomes  known,  the  scheduler  can  transfer 
those  aircraft  from  routes  where  demand  is  less  than  anticipated  to 
those  where  it  is  higher  than  expected.  Thus,  the  scheduling  problem 
is  a  dynamic  one  of  periodically  reassigning  capacity  to  routes  in 
order  to  minimize  the  expected  costs  of  the  operation  over  time. 

Assigning  insufficient  airlift  increase?  both  the  costs  due  to  delay 
and  the  procurement  costs  of  "emergency"  commercial  lift;  assigning 
excess  capacity  results  in  an  opportunity  cost  of  foregoing  service 
elsewhere  in  the  system.  Since  both  situations  are  stochastic,  the 
best  that  can  be  done  is  to  minimize  the  cost  entailed  in  these  situ¬ 
ations  on  the  average. 

BASIC  APPROACH 

We  turn  now  to  an  overview  of  our  strategy  in  treating  the  sched¬ 
uling  problem.  Ideally,  the  problem  should  be  treated  as  a  multiple 
period  (30  days  or  longer),  sequential  decision  problem  with  account 
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being  taken  of  daily  demand  variation  and  physical  constraints  on  air¬ 
craft  in  each  of  the  ensuing  periods  in  the  current  horizon.  Such  an 
approach  is  clearly  not  feasible  fr< a  a  computational  standpoint  for  a 
system  the  size  of  MAC  or  even  each  separate  Air  Force.  Therefore, 
broadly  following  the  current  procedures  at  MAC,  we  break  the  process 
into  two  parts  --  monthly  and  daily. 

A  monthly  schedule  is  determined  that  accounts  for  uncertainty  for 
the  forthcoming  month  considered  as  one  aggregate  period  in  which  deci¬ 
sions  planned  in  advance  may  be  modified.  Although  this  is  far  from 
ideal,  it  is  superior  to  ignoring  uncertainty  and  the  subsequent  oppor¬ 
tunity  to  remedy  its  effects.  The  second-stage  decisions,  as  we  shall 
see,  are  "virtual"  or  fictitious,  but  are  nonetheless  important  because 
they  help  shape  a  good  first-stage  decision,  the  initial  assignments 
that  serve  as  the  starting  point  for  the  daily  modifications. 

The  monthly  flights  are  apportioned  to  the  days.  For  example,  if 
fifteen  flights  of  a  given  aircraft  are  found  from  the  monthly  model, 
this  may  be  interpreted  as  one  flight  per  aircraft  every  two  days  during 
the  month  if  the  given  aircraft  type  can  fly  the  mission. 

Once  the  month  commences,  daily  amendments  to  the  number  of  flights 
scheduled  for  that  particular  day  are  determined  by  the  daily  model, 
which  is  also  a  two-stage  optimization.  Changes  in  assignments  optimal 
for  the  forthcoming  two-day  period  are  determined.  They  minimize  the 
sum  of  costs  incurred  on  the  current  day,  given  the  amount  of  cargo  on 
hand  plus  the  expected  value  of  tomorrow's  costs,  which  depend  upon 
both  the  amount  of  new  cargo  and  the  changes  made  on  the  first  day. 
Consecutive  decisions  of  the  daily  model  are  thereby  chained  together 
in  a  recursive  fashion  with  outputs  from  the  model  for  day  t  becoming 
inputs  to  the  model  for  day  t+1. 
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III.  MATHEMATICAL  FORMULATION 


In  this  section  we  formulate  the  monthly  and  daily  models  as  two- 
stage  linear  programs  under  uncertainty  [2],  [9]  with  a  random  right- 
hand  side  in  the  second  stage.  Several  of  the  program's  variables 
represent  numbers  of  different  kinds  of  flights  and  consequently  must 
be  integer  valued.  However,  this  aspect  of  the  problem  is  neglected 
in  order  to  obtain  a  program  that  is  computationally  solvable,  thus 
leaving  the  user  to  adopt  a  rounding  scheme.  Ihe  general  form  for 
problems  of  this  type  is  the  following: 

^  min  -jcx  +  E  [min  qy  |Wy  =  E  -  Tx]T  subject  to  Ax  *  b. 
xaO  1  »  ySO  J 

In  this  problem  c,  q,  and  b  are  known  vectors  of  dimension  n^ ,  n^,  and 
m^  respectively.  A,  T,  and  W  are  known  matrices  of  dimension  m^  x  n^ , 
m2  x  n^ ,  and  m^  x  n2»  respectively.  The  vectors  x  and  y  represent  the 
first-  and  second-stage  decision  vectors.  §  is  a  random  vector  with  a 
known  probability  distribution  and  E(*)  denotes  the  expectation  with 
respect  to  that  distribution. 

The  problem  can  be  given  the  following  interpretation.  Before 

the  random  variable  £  is  observed,  one  must  choose  nonnegative  values 

of  x  satisfying  Ax  «  b.  The  immediate  cost  of  this  choice  is  cx.  After 

x  is  chosen  and  the  random  variable  t  is  observed,  nonnegative  values 

must  be  chosen  for  y  which  satisfy  Wy  =  f  -  Tx  at  a  cost  of  qy.  The 

chosen  value  of  y  will,  for  ecch  observed  value  of  §,  be  the  one  which 

minimizes  qy.  Letting  Q(x)  =  E  [rain  qy  |Wy  =  §  -  Tx},  we  wish  to  minimize 

5  ysO 

the  first-stage  cost,  cx,  plus  the  expected  second-stage  cost,  Q(x) . 

The  problem  can  be  viewed  as  a  dynamic  programming  problem  with  only 
two  stages,  with  y  a  function  of  the  state  vector  5  “  Tx.  Note  that 
x  does  not  depend  on  the  observed  values  of  §  but  only  on  the  distribu¬ 
tion  of  §,  while  y  does  depend  on  the  observed  values  of  §  and  on  the 
chosen  values  of  x. 

MONTHLY  MODEL 


In  the  monthly  model  the  sequence  of  decisions  is  the  following: 
1.  The  number  of  flights  of  each  aircraft  type  in  the  MAC  fleet 


-6- 


*  Vert 

over  each  route  is  assigned  before  monthly  requirements  on 
each  route  are  known  with  certainty. 

2.  After  monthly  -equirements  are  observed,  some  flights  assigned 
to  routes  with  lower  than  expected  demand  are  switched  to 
routes  with  higher  than  expected  demand.  Commercial  airlift 
is  also  sometimes  added  to  routes  with  excess  demand. 

The  monthly  model  has  the  following  structure.  The  first-stage 
constraints  state  that  for  each  aircraft  type  the  total  number  of  flyi.ng 
hours  allocated  to  all  routes  cannot  exceed  the  total  number  of  flying 
hours  available  of  that  type.  The  second-stage  constraints  are  of  two 
types.  The  first  specifies  that  the  number  of  flying  hours  of  a  given 
aircraft  type  diverted  from  a  particular  route  to  other  routes  cannot 
exceed  those  initially  assigned  to  it.  The  second  type  are  demand 
balance  equations,  which  state  that  for  each  route  the  total  carrying 
capacity  (that  originally  assigned,  minus  that  diverted  to  other  routes 
plus  that  diverted  from  other  routes)  minus  unused  carrying  capacity 
plus  unsatisfied  demand  is  equal  to  total  demand  for  that  route. 

The  objective  function  consists  of  the  cost  of  the  final  flying 
program  (the  initial  plus  the  amended  assignments)  plus  penalty  costs 
of  excess  demand  or  supply.  The  cost  of  excess  demand  is  reflected  in 
both  the  cost  of  additional  commercial  lift  plus  the  extra  flying  time 
consumed  in  switching  aircraft  from  one  route  to  another.  Specifically, 
the  program  is  as  follows. 


Find  min  Z,  x. . ,x. , y+,y. 

’  ij  ijk’-'j’-'j 


such  that 


A  route  is  a  channel  or  SAAM  mission  from  each  AFOE.  Routes  are 
assumed  to  be  point-to-point  with  no  intermediate  on-load  or  off-load 
of  cargo.  However,  e.g.,  a  trip  from  Travis  to  Bangkok  via  Hickam  is 
considered  a  different  route  from  a  flight  which  goes  via  Guam. 
rkic 

Requirements  as  used  here  will  mean  net  requirements  from  which 
capacity  of  commercial  airlift  on  an  annual  contract  (for  which  the 
frequency  of  trips  and  routes  is  fixed)  and  any  other  fixed  capability, 
such  as  Air  National  Guard  flights,  is  subtracted.  This  has  the  effect 
of  merely  lowering  the  mean  of  the  distribution  of  requirements. 
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<2)  IcLjxij +  E{min[  I  (ct.k  -  xijk  +  IcM +  Z cjyJ]} =  z 

i,j  i-J.kyj 


tj 


First  stage: 
(2.1) 


/  a . .x. .  £  F. ,  alii 
Li.  j  tj  l’ 


J 


Second  stage: 

(2.2)  x.  .  -  Y  x. . .  2t  0 ,  all  i,j 

ti  L.  a..  iik 


k^j 


LJ 


(2.3)  Y  b  x.  .  -  y  k.  ( x...+  y  b.  .x. .  .  +  y*  -  y .  =  d.  ,  al  1  j  , 

o  ij  i]  L  L  ij '  a. ./ J  Ljk  L  i]  Lkj  j’ 

i  L,k/j  J  i,kjij 


where  x. . 

1J 


x.  .. 
ijk 


+ 

yj 


and 


number  of  flights  of  aircraft  type  i  initially  assigned 
to  route  j  during  the  month; 

number  of  flights  of  aircraft  type  i  assigned  to  route  k 
using  hours  made  available  by  canceling  route  j  flights; 

demand  on  route  j  which  is  satisfied  by  commercial  lift, 
if  permitted,  or  unsatisfied  demand  if  commercial  lift 
not  permitted; 

unused  capacity  on  route  j; 


F. 

i 


d. 

J 


number  of  hours  required  by  aircraft  type  i  for  a  flight 
initially  assigned  to  and  flown  on  route  j; 

number  of  hours  required  by  aircraft  type  i  for  a  flight 
on  route  k  that  uses  hours  made  available  by  canceling 
route  j  flights  (aijk  a  a.k)  ; 

the  carrying  capacity  (tons  or  any  other  appropriate 
measure'  of  a  flight  of  an  aircraft  of  type  i  on  route  j; 

the  maximum  number  of  flying  hours  for  aircraft  of  type  i 
available  during  the  month; 

total  demand  (tons  or  any  other  measure)  for  route  j  (a 
random  variable); 


t 
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ij 


cijk 


cost  of  flight  of  aircraft  type  i  initially  assigned  to 
and  flown  on  route  j  (c_  s  0); 

cost  per  flight  of  aircraft  type  i  assigned  to  route  k  from 


hours  made  available  by  canceling  route  j  flights  (c..,  >  c,,) 

1 J  K  IK 


cost  per  ton  of  commercial  augmentation  on  route  j.  If 
commercial  augmentation  is  not  available  to  carry  excess 
demand,  this  may  instead  represent  a  shortage  cost; 


J 


cost  of  a  unit  of  unused  carrying  capacity  on  route  j. 


Note  that  since  a  flLght  assigned  to  and  flown  on  route  j  takes 
a  hours  while  a  flight  flown  on  route  k  from  hours  diverted  from 


route  j  takes  a£jjc  hours,  such  a  flight  on  route  k  results  in  the  can¬ 


cellation  of  ^a£j^/a£j)  flights  on  route  j.  Thus,  of  the  flights  ini¬ 


tially  assigned  to  route  j,  those  th  t  are  actually  flown  number 


x. 


k*j 


ijk 


This  is  reflected  in  expressions  (2),  (2.2),  and  (2.3). 

The  first-stage  constraints  correspond  to  Ax  =  b  and  the  second- 

stage  constraints  to  Wy  *  5  "  Tx  in  expression  (1).  Notice  that  the 

constraints  (2.2)  are  placed  in  the  second  stage  despite  the  fact  that 

the  right-hand  side  is  not  a  random  vector.  This  is  done  because  the 

switching  variables  x..,  are  assumed  to  be  determined  after  demand 

L  j  K 

becomes  known  and  thus  they  depend  upon  information  about  the  d  .  Also, 

+  J 

they  are  to  be  chosen  simultaneously  with  the  y^  and  y^  to  minimize 
the  overall  second-stage  cost  function. 


DAILY  SCHEDULING  MODEL 

The  daily  scheduling  model  is  also  a  two-stage  linear  program 
under  uncertainty.  Its  principal  first-stage  decision  variables  are 
the  number  of  flights  of  each  type  to  switch  from  one  route  to  another 
in  the  forthcoming  day  —  variables  that  played  the  role  of  second- 
stage  recourse  variables  in  the  monthly  model.  In  addition,  the  daily 
model  is  couched  in  terras  of  unmoved  cargo  stocked  at  the  various  origins 
(aerial  ports  of  embarkation)  for  various  lengths  of  time;  the  <_argo  has 
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not  yet  been  shipped,  and  subsequently  incurs  certain  delay  costs. 
Consequently,  the  amount  of  cargo  in  storage  for  various  lengths  of 
time  waiting  to  be  loaded  on  different  aircraft  types  over  the  differ¬ 
ent  routes  constitutes  another  set  of  first-stage  decision  variables. 

In  the  first  stage  the  constraints  consist  of  balance  limitations 
on  the  number  of  flying  hours,  bounds  on  switches,  capacity  restrictions 
on  shipments,  and  inventory  balance  equations  for  cargo  which  has  been 
in  the  port  for  over  one  day  —  quantities  known  with  certainty.  The 
second  stage  consists  of  demand  equations  for  newly  arriving  cargo  which 
is  treated  as  a  random  variable.  These  equations  include  terms  to 
account  for  a  possible  under-  or  over-supply  of  capacity.  It  will  be 
assumed  that  for  a  given  route  new  cargo  has  a  lower  priority  (i.e.  , 
lower  penalty  or  delay  cost)  than  old  cargo.  The  reason  for  this  will 
become  clear  later. 

The  objective  function  is  again  one  of  cost  minimization.  In  this 
case  the  costs  include  (1)  incremental  operating  costs  resulting  from 
flight  changes,  (2)  daily  commercial  procurement,  (3)  delay  costs  for 
cargo  that  remains  unmoved  for  various  lengths  of  time,  and  (4)  expected 
penalty  costs  resulting  from  misallocation  of  capacity. 

To  write  the  daily  model,  we  introduce  the  following  notation 
(variables  ref?r  to  the  t  day  of  the  month): 

x  “  number  of  flights  of  aircraft  type  i  previously  scheduled 
J  to  fly  over  route  j,  a  known  constant; 

X.,,  *  number  of  flights  of  aircraft  type  i  switched  to  route  k 

J  using  hours  made  available  by  canceling  route  j  flights; 

yf  *  capacity  reserved  for  cargo  £  days  old  shipped  on  aircraft 
^  type  i  over  route  j  (£=0,1,...,L); 

w^.  =  capacity  reserved  for  cargo  £  days  old  shipped  on  commer- 
J  cial  lift  of  type  h  over  route  j; 

z.  t  *  number  of  flights  of  commercial  lift  of  type  h  procured  for 
“J  route  j; 

s^  *  amount  of  cargo  i  days  old  remaining  to  be  shipped  over 
■*  route  j  at  the  end  of  the  day; 

l 

H  *  inventory  of  cargo  £  days  old  for  route  j  for  shipment  at 

J  the  beginning  of  the  day  (H°  represents  new  cargo  and  is 

a  random  variable) ;  ^ 

v^  *  unused  capacity  on  route  j  at  the  end  of  the  day; 


i  i 
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=  number  of  hours  required  by  aircraft  type  i  scheduled  and 
flown  on  route  j; 

=  number  of  iiours  required  by  aircraft  type  i  for  a  flight 
on  route  k  using  hours  made  available  by  canceling  route 
J  flights  (atJki  aik); 

«*  the  carrying  capacity  of  aircraft  type  i  over  route  j; 

=  cost  per  flight  of  aircraft  type  i  scheduled  and  flown  on 

route  j  (c .  .  «*  0) : 
ij 

=  cost  per  flight  of  aircraft  type  i  assigned  to  route  k  using 
hours  made  available  by  canceling  route  j  flights  (c  ^  £  c  )  ; 

=  cost  per  flight  of  commercial  lift  of  type  h  over  route  j; 

=  unit  cost  per  ton  of  one  additional  day's  delay  for  cargo 
l  days  old  which  is  to  be  shipped  over  route  j.  It  is 
assumed  that  v°  =  min 

J  l  J 

=  unused  capacity  on  route  j. 


Then  the  problem  for  day  t  can  be  stated.  Find 


Mia  Z.xljk..J,*hJ.yfJ.»Jj.vJ  *0 


such  that 


i  (■ 

i,j  ,k^j  \ 


uk  -  'ijC^Kjk +  I 

ij  /  j 


+  Y  8.  .  z.  .  +  E-jmin  /  y°,s°.  +  )  q.v.l  3  Z> 

L  hj  hj  l  „  Uy j  j  j  jJ 


Uj’Vj  j 


First  Stage: 
(3.1) 


Y  — Liii  x. s  x. . , 
a.,  «-JK  LJ 

Mj  tj 


all  i  and  j 


(3.2)  Xyfj  -  yij  +  I  (bij^ikjk  -  Ibii*ikJ  5  Vis'  311  ‘ and  1 

I'O  k<j  aij  Mj 
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(3.3)  IuSj  +“hj  -  -hj=hj  5  °’  aUhandj 

IX 

(3«4)  Y  yf.  +  )  wf .  +  *  H^,  alljand£>0- 

4-  ij  w  h]  j  J 

i  h 


Second  Stage: 


(3.5) 


all  j . 


Equation  (3.1)  states,  as  in  the  monthly  model,  that  the  number 
of  flying  hours  switched  off  a  given  route  and  aircraft  type  cannot 
exceed  those  originally  assigned. 

The  meaning  of  Eqs.  (3-2)  and  (3-3)  is  that  the  total  tonnage 
loaded  of  cargo  already  received  plus  the  amount  reserved  for  new  cargo 
on  a  given  aircraft  type  over  a  given  route  must  be  less  than  the  carry¬ 
ing  capacity  of  all  aircraft  of  chat  type  assigned  to  the  route  for  MAC 
and  commercial  vehicles,  respectively.  Equations  (3-4)  state  for  each 
duration  of  cargo  destined  for  each  route,  the  total  tonnage  allocated 
for  MAC  aircraft  and  commercial  aircraft  plus  the  amount  not  moved  must 
equal  the  inventory  of  cargo  at  the  beginning  of  the  day  (a  known  con¬ 
stant  for  all  but  newly  arriving  cargo).  We  note  that  the  assumption 
that  newly  arrived  cargo  has  lowest  penalty  cost  (mathematically  stated, 

y°  -  min  (yf)  for  all  j)  assures  that  no  new  cargo  will  be  carried  unless 
J  £  J 


all  old  cargo  is  carried.  This  assures  that  for  all  but  newly  arrived 
cargo,  the  amount  shipped  uses  all  Che  allotted  capacity.  If  this 
assumption  were  not  made,  it  would  be  more  reasonable  to  determine  the 
y°^  and  w?^  simultaneously  with  the  >  0)  as  a  function  of 

the  actual  amount  ol  new  cargo.  This  would  necessitate  shifting  equa¬ 
tions  (3>2)  through  (3*4)  into  the  second  stage,  making  for  a  large 
linear  program  for  the  second  stage.  The  present  formulation  permits 


★ 

For  routes  on  which  commercial  augmentation  is  not  permitted, 
the  corresponding  z^.  are  omitted. 
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t'ne  second  stage  to  have  a  very  simple  structure  whLch  reduces  the  com¬ 
putational  burden  in  solving  it.  This  matter  wiii  be  discussed  In  more 
detail  below. 

Finally,  Eq.  (3.5)  states  that  for  each  route  the  inflow  of  new 
cargo  minus  the  total  capacity  reserved  for  shipment  of  new  cargo  on 
all  types  of  MAC  and  commercial  vehicles  is  equal  to  any  cargo  not 
shipped  minus  any  unused  capacity. 

In  addition,  we  may  wish  to  include  constraints  that  do  not  satu¬ 
rate  the  capacity  of  enroute  bases.  Let  J(p)  equal  the  set  of  routes 
which  use  base  p  as  an  enroute  stop  for  refueling  or  crew  change.  For 
example,  channels  from  Norton  to  Saigon  via  Hickam  would  be 
routes  with  Hickam  as  base  p.  Let  ^  be  the  maximum  number  of  flights 
that  can  be  accommodated  at  base  p  in  a  single  day.  Then  these  con¬ 
straints  can  be  expressed  as 

I  (xtj  -  I  “jfr  *ijk  •  1  *Lkj) 1  °P  • 

i»jeJ(p)  k^j  2  k^j 


l 

Letting  Hj(t)  be  the  inventory  of  cargo  for  route  j  that  is  i  days 

old  on  day  t,  and  letting  quantities  in  parenthesis  show  time  dependence 

l  l- 1 

of  other  quantities  in  a  similar  way,  note  that  H,(t  +  1)  =  s  (t)  ,  using 

£-1  2  2 
the  convention  that  s  (0)  refers  to  cargo  not  shipped  on  the  last  day  of 

the  preceding  month. 


In  addition. 


the  quantities  x 


ij 


are  updated  daily  by  the  relations 


where  the  Xj^  and  are  values  of  the  variables  as  deter¬ 

mined  by  the  model  for  day  t. 

In  the  monthly  model  there  are  m  (number  of  Individual  aircraft 


types)  constraints  in  the  first  stage  and  (mn  +  n)  in  the  second  stage 
(n  being  the  total  numoer  of  routes).  Assuming  five  aircraft  types* 


* 


Currently  MAC  uses  the  C-124,  C-130,  C-133,  C-135  and  C-141. 
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and  100  routes,  this  presents  a  linear  program  for  the  first  stage  with 
five  constraints  and  600  for  the  second  stage.  However,  in  the  second 
stage  the  first  mn  are  of  the  generalized  uppefbound  type,  to  be  dis¬ 
cussed  later,  to  which  special  methods  may  be  applied. 

In  the  daily  model  there  are  at  most  n(2m  +  H  +  L)  constraints  in 
the  first  stage  and  n  itt  the  second  stage,  where  H  is  the  number  of  indi¬ 
vidual  commercial  vehicles  and  L  is  the  number  of  duration  classes  of 
old  cargo.  For  H“3,L=3,m=5,  n  =  100,  this  yields  one  first- 
stage  program  with  1600  constraints.  The  second  stage  contains  100 
equations.  While  this  is  within  the  range  of  linear  programming  codes 
for  "tnird  generation"  computers,  the  first-stage  problem  can  be  sim¬ 
plified  by  treating  if  by  decomposition  methods  where  Eq.  (3.2)  are 
the  coupling  equations  and  (3.1),  (3. 3-3. 4)  are  the  two  disjoint  sub¬ 
problems.  Again,  for  the  former  subproblem,  the  majority  of  the  con¬ 
straints  (3.1)  are  of  the  generalized  upper-bound  variety. 


miipi 


vsammBNmm  .j ,.jpn.»  im  i 
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IV.  SOLUTION  TECHNIQUES 

This  section  describes  "wo  alternative  computational  methods  for 
solving  two-stage  linear  programs  under  uncertainty.  These  methods  may 
be  applied  tc  the  monthly  and  daily  model.  They  assume  the  random  vari¬ 
ables  have  a  continuous  distribution.  Each  is  a  variant  of  an  algorithm 
for  convex  programming,  which  is  required  for  this  type,  of  problem.  The 
first  is  a  cutting-plane  or  constraint-generating  procedure,  while  the 
second  is  a  gradient  method.  We  outline  the  steps  required  for  the 
general  problem  exemplified  by  the  monthly  model  and  indicate  certain 
simplifications  that  can  be  made  in  the  daily  model  whose  second-stage 
has  a  special  str~_ture.  Before  proceeding  to  the  algorithms,  we  review 
some  salient  facts  concerning  two-stage  linear  programming  under  uncer¬ 
tainty. 

BASIC  FACTS  FROM  LINEAR  PROGRAMMING  UNDER  UNCERTAINTY 

As  mentioned  in  Sec.  Ill,  both  the  monthly  and  daily  problems  are 
special  cases  of  the  problem: 

(4)  Find  x  2:  0,  min  Z  such  that 

Z  =  cx  +  0(x) 

Ax  =  b 

where  Q(x)  =  E  q>  ,  'y  =  §  -  TxJ  • 

Wets  [9]  has  shown  that  Q(x)  is  convex  in  x.  Thus  the  above  linear 
program  under  uncertainty  i?  equivalent  to  a  deterministic  convex 
prog"  m  in  x.  The  algorithms  to  be  developed  focus  and  are  based  on 
certain  properties  of  the  function  Q(x). 

One  important  property  (to  be  proved  later  in  the  Constraint  Gen¬ 
eration  section)  is  that  the  linear  function  p(x°)  -  (tt(x0)T)x0  is  a 

support  function'  (or  tangent  plane)  to  Q(x)  at  x°eK,  where  n(x°)  is  the 
expectation  of  the  optimal  dual  variables  and  p(x°)  is  the  expectation 


A  linear  function  li(x)  is  a  support  for  convex  function  z(x)  de 
fined  on  XcEn  at  x°  if  z(x)  i  H(x)  for  all  xeX  and  z(x°)  =  H(x°) . 
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of  the  demand  terms  Ln  the  dual  objective  function.  Specifically,  let 
rt(x°,5)  be  an  optimal  solution  to 

(5)  -  Tx°) |rfr  S  q}  . 

Then  n(x°)  =  Ec.rr(x°,f)  and  p(x°)  =  E^fn(x°,§)|]  , 

This  property  forms  the  basis  for  the  cutting-plane  method. 

Another  important  property  [9]  is  that  if  the  distribution  of  c 
is  absolutely  continuous,  Q(x)  is  differentiable  at  x°eK  and  has  the 
gradient  -rr(x°)T,  where  again  tt(x°)  is  the  expected  value  of  the  solu¬ 
tion  to  (5).  This  property  is  utilized  in  the  gradient  approach  de¬ 
veloped  by  Wets  [8]. 

In  general,  the  -f(x°)  and  p(x°)  cannot  be  explicitly  computed; 
consequently,  the  methods  to  be  described  resort  to  sampling  procedures 
to  estimate  these  quantities. 

Our  monthly  and  daily  problems  both  possess  the  property  of 
complete  recourse.  In  other  words,  each  x  satisfying  Ax  =  b  yields  a 
feasible  program  for  all  possible  E  in  the  second-stage  constraints 
Wy  =  5  *  Tx.  This  permits  us  to  eliminate  one  of  the  steps  in  the 
algorithms  to  be  described  for  finding  feasibility  cuts  [7,  pp.  31-34]. 
Moreover,  che  daily  model  exhibits  simple  recourse,  i.e.,  W  =  [i,  -ij. 
This  permits  simplification  of  one  major  step  in  the  following  algorithm. 

CUTTING-PLANE  METHOD 

A  cutting-plane  algorithm  has  been  developed  by  Wets  and  Van  Slyke 
[7]  for  solving  the  general  problem.  It  can  be  viewed  as  an  application 
of  Kelley’s  cutting-plane  method  [5]  to  the  equivalent  convex  program. 

In  this  approach  the  convex  function  Q(x)  appearing  in  the  objec¬ 
tive  function  is  treated  as  a  constraint.  The  objective  function  is 
replaced  by  z(x)  =  cx  +  9  and  the  constraint. 


(6) 


Q(x)  s  9, 
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considere  the  first-stage  program  where  6  is  treated  as  an  addi¬ 
tional  decision  variable.  Since  the  second-stage  costs  are  nonnegative 
in  our  problem,  we  may  consider  9  as  a  nonnegative  variable. 

The  constraint  (6)  is  represented  by  a  set  of  support  planes 
1c  k 

p  -  (tt  T)x;  i.e.,  a  convex  function  is  represented  as  tl  upper  enve¬ 
lope  of  its  supporting  planes.  Four  constraints  are  generated  by  solv¬ 
ing  the  second-stage  problem  for  a  fixed  x  which  solves  an  augmented 
first-stage  problem.  Thus  the  technique  consists  of  a  decomposition 
of  the  overall  problem  into  two  linear-programming  subproblems.  The 
method  can  be  viewed,  in  fact,  as  a  particular  case  of  Benders'  par¬ 
titioning  method  [l]  applied  to  the  original  problem. 

The  particular  problem  that  is  solved  at  any  iteration,  say  the 
th 

k  ,  is  the  following: 

1.  Find  x,9  2  0,  min  z"  such  that 

(7)  cx  +  9  -2 

a)  Ax  *  b 

b)  (rrkT)x  +  0  *  pk,  k  =  1 . S, 

where  S  is  an  index  of  a  subset  of  previously  generated  constraints  to 
be  described  more  fully  below. 

The  constraints  (7b)  are  generated  sequentially  and  have  the  pro¬ 
perty  that  for  all  x  satisfying  (7a),  (x,Q(x))  =  (x,8)  satisfies  (7b). 

2.  Let  (x,G)  be  the  solution  obtained  from  step  l  and  solve  for 
Q(x) .  If  Q(x)  =  6,  x  solves  (A).  Otherwise,  a  constraint  of  the  type 
(7b)  is  generated  that  is  violated  by  (x,9). 

Constraint  Generation 

For  fixed  x  =  x  and  a  given  §  the  second-stage  linear  program 

becomes 

(8)  Find  y  2  0,  0(x,§)  where 

Q(x,D  =  min  qy  subject  to  Wy  =  ?  -  Tx. 
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Let  n(x,§)  be  the  simplex  multipliers  corresponding  to  an  optimal 
basis  in  (8).  It  follows  that  n(x,§)(§  -  Tx)  =  Q(x,5)  and  tt(x,^)W  £  q. 
Thus,  n(x,§)  is  feasible  in  the  following  program  for  all  x: 

(9)  Find  tt.  max  w  where 

n(5  -  Tx)  =  w 
TtW  s  q  . 

Since  (9)  is  dual  to  (8)  for  general  x,  it  follows  from  duality 
theory  that 

n(*,l)(C  -  Tx)  *  Q(k,S) 

for  all  x.  Letting  tt(x)  =  E-rr(x,t)  and  p(x)  =  EpC^x , §)  E] ,  it  follows 
that  for  all  x,  (n(x)T)x  +  Q(x)  2  p  (x)  .  Hence,  the  following  constraint 
must  be  satisfied  by  all  (x,Q(x)): 

(10)  (n(x)T)x  +  0  a  o(x). 

Algorithm 

All  constraints  of  the  type  (7b)  are  of  type  (10).  That  is,  each 
k  k 

(n  ,P  )  is  equal  to  some  (tt(x)  , P  (x) )  .  Since  for  every  x  satisfying 
Ax  *  b,  (x,Q(x))  is  feasible  for  (7),  it  follows  that  in  (7),  min  z  s 
min  z,  where  min  z  solves  (4).  Thus,  if  (x,0)  solves  (7)  and  Q(x'  -  9, 
cx  +  Q(X)  £  min  z  and  X  solves  (4)  . 

On  the  other  hand,  if  Q(x)  f-  0,  it  follows  (since  cx  +  Q(x)  ^ 
min  z)  that  Q(x)  >  0.  Also  (rr(x)T)x  +  Q(x)  =  p(x)  by  taking  expecta¬ 
tions  in  the  expression  tt(x,§)(§  -  Tx)  =  Q(x,£).  Thus  (x,0)  violates 
the  expression  (10) . 

One  problem  which  arises  is  that  it  is  possible  to  generate  a  large 
number  of  constraints  of  type  (7b).  However,  in  reference  [6]  it  is 
shown  that  all  constraints  in  (7b)  that  are  not  tight  (i.e.,  where 
equality  does  not  hold)  in  the  optimal  solution  to  (7)  may  be  dropped 
from  succeeding  iterations.  Thus,  the  number  of  such  constraints  never 
exceeds  the  number  of  rows  in  T.  Formally,  the  algorithm  is  as  follows 
with  S  initially  zero  in  step  l. 


I 
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1.  Solve  for  x  i  0,  min  z  such  that 
cx  +  e  =  ~z 

Ax  =  b 

(nkT)x  +  0  -  sk  =  pk  k  -  1,...  ,s. 

Let  x  be  the  optimal  x  obtained. 

2.  If  s  >0,  delete  the  constraint  in  which  sk  appears,  k  =  1, 

• • . >S  and  renumber  those  remaining. 

3.  Sample  from  the  distribution  of  %.  For  each  point  gl(i  = 

...  ,n)  in  the  sample,  solve  for 

y  ^  0,Q(x,D  where 
Q(x,?L)  =  min  qy  subject  to 
Wy  =  ?L  -  1*. 

Increase  S  by  1  and  let  ttS  =  n'^nCx,?1)  and  pS  =  n'lE?lTr(x , £l)  ,  where 

n  is  the  sample  size  and  TT(x,§i)  are  the  optimal  simplex  multipliers 
for  the  sample  point 

4.  If  Q(x)  -  0  s  6,  where  6  is  a  predetermined  constant  or  vari¬ 
able,  terminate.  Otherwise,  return  to  l. 

The  quantity  6  in  step  4  may  be  chosen  by  any  sensible  criteria. 

For  example,  it  may  be  a  constant  or  it  could  be  a  small  fraction  times 
the  last  step  1  solution  (i.e.,  z(x)). 

When  the  second  stage  exhibits  simple  recourse  (i.e.  ,  the  only 
second-stage  variables  are  slack  variables)  as  is  *  .e  case  in  the  daily 
model,  step  3  can  be  considerably  simplified.  In  this  case  the  sampling 
procedure  and  the  need  to  solve  a  linear  program  can  be  dispensed  and  the 
expected  value  of  the  shadow  prices  calculated  explicitly.  Specifically, 
Wy  =  -  Tx  may  be  expressed  as  q^  -  q"  =  e1  -  Tx,  i  =  l,...,m  , 

and  from  [8] 
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That  is,  they  are  equal  to  the  cost  of  the  iC^  excess  demand  times  the 
probability  of  its  occurrence  minus  the  cost  of  the  excess  supply  in 
the  it}l  equation  times  the  probability  of  its  occurrence. 

The  Second-Stage  Problem 

The  sampling  involved  in  the  calculation  of  Q(x)  in  step  3  of  the 
algorithm  involves  the  solution  of  many  linear  programs.  The  solution 
time  for  this  step  can  greatly  be  reduced  by  taking  advantage  of  the 
fact  that  if  a  basis  is  optimal  for  a  program  and  the  right-hand  side 
is  changed,  that  basis  will  either  be  optimal  or  infeasible  for  the 
new  right-hand  side.  Thus,  the  sequence  of  programs  for  the  second 
stage  may  be  sc’ved  as  follows.  Solve  the  program  for  one  sample  point. 
Check  each  of  the  other  sample  points  for  feasibility  with  respect  to 
the  optimal  basis  found.  (The  rr  vector  associated  with  this  basis  is 
an  optimal  one  for  all  sample  points  that  are  feasible.)  Then  solve 
the  program  for  some  sample  point  that  was  infeasible  and  repeat.  This 
procedure  reduces  the  number  of  linear  programs  from  the  number  of  ele¬ 
ments  in  the  sample  to  the  number  of  optimal  bases  for  the  sample  points. 
After  one  optimal  basis  is  found,  others  may  of  course  be  found  from  this 
one  by  a  sequence  of  dual  simplex  pivot  steps.  Hence,  each  pivot  will 
yield  a  basis  whose  cost  *-ow  Is  nonnegative.  It  is  suggested,  then,  that 
after  each  iteration,  fe.  bility  be  checked  for  the  other  sample  points 
which  have  not  been  solved  or,  rather  than  wait  until  the  current  one 
being  considered  is  feasible.  If  is  also  suggested  that  the  optimal 
basis  be  saved  from  iteration  to  i.'eration. 

GRADIENT  METHOD 

In  this  approach  the  equivalent  convex  program  is  solved  by  a  non¬ 
linear  programming  algorithm  which  utilizes  the  gradient  of  the  objective 
function  of  (A)  at  x° , 

y  z(x°)  =  (c  -  tt(x°)T)  . 

The  method  we  describe  is  based  upon  a  variant  of  the  Frank-Wolfe  algo¬ 
rithm  [4]  for  convex  programming  with  linear  constraints  developed  by 
Wets  [8]  for  stochastic  programs.  Other  gradient  methods  for  convex 
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programming  such  as  the  Gradient  Projection  method,  Zangwiil's  convex 
simplex  method  [10]  or  Zoutendijk1 a  method  of  feasible  directions  [11] 
could  be  employed  as  well. 

In  brief,  this  method  obtains  solutions  to  a  sequence  of  linear 
approximations  to  the  given  problem  of  the  form 


(ID 


min  Vc(x^  ^)x 
x 


subject  to 

Ax  =  b 
x  2  0  , 


£_  1 

where  x  is  the  solution  for  x  to  the  given  problem  at  the  preceding 

£ 

or  (£-l)st  iteration.  An  improved  feasible  solution  x  is  obtained  by 
£_!  £  £_1  _£ 
adding  to  x  an  improving  feasible  direction  \(x  -  x  ),  where  x 

is  the  (basic)  optimal  solution  to  (11)  and  \  (a  scalar)  is  determined 

by  solving  a  one-dimensional  optimization  problem  to  be  described  below. 

Several  features  of  this  method  should  be  noted.  First,  the  linear 
program  which  is  actually  solved  in  the  first  stage  is  not  augmented 
but  contains  only  constraints  of  the  first  stage;  only  the  cost  func¬ 
tion  varies  from  iteration  to  iteration.  Second,  the  computational 
labor  is  again  split  between  first  and  second  stages  with  the  expected 
values  of  the  shadow  prices,  which  are  used  in  the  objective  function 
of  the  first-stage  problem,  generated  in  a  sampling  procedure  in  solving 
a  set  of  second-stage  programs  for  a  fixed  x.  Third,  since  the  method 
is  a  gradient  approach,  it  is  self-correcting  in  the  sense  that  if  a 
poor  approximation  to  the  gradient  is  generated  which  leads  to  a  small 
improvement  in  the  solution,  it  is  likely  to  be  corrected  at  later  itera¬ 
tions.  By  contrast,  in  the  cutting-plane  approach  it  is  possible  that 
sampling  variations  occasionally  generate  support  functions  which  cut 
inside  the  feasible  region  and  indicate  falsely  that  an  optimum  has  been 
reached. 
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Aleorithm 

The  exact  steps  of  the  algorithm  are  as  follows: 

1.  (Initialization.)  Set  iteration  count  £  =  0  and  tt  =  ff(x  ) 

=  0  [the  sample  estimate  of  the  expected  value  of  the  optimal  multi¬ 
pliers  for  the  second-stage  program]. 

2.  (Generate  optimal  solution  to  the  linearized  equivalent  convex 
program.)  Set  £  =  £  +  1.  For  £  =  1,2,...,  solve  the  linear  program 


(0 


min  (c  -  fr^  ^T)x 
x 


subject  to 


Ax  =  b 
x  S  0. 


a 

Denote  the  optimal  solution  to  (i)  by  x  .  For  £  =  1 ,  go  to  step  5. 
Otherwise  go  to  step  3. 

3.  (Seek  improved  direction.)  Solve  the  scalar  optimization 
problem 


(ii)  min  \jr(X)  =  c[(l  -  \)xl  1  +  \x^] 

OA^i 

+  Q[(i  -  \)xl  1  +  Xx*], 


where  Q(x)  is  a  sample  estimate  of  Q(x)  obtained  by  solving  the  second- 
stage  program  for  a  random  sample  of  l  and  a  fixed  x.  In  (ii)  this  func¬ 
tion  is  viewed  as  a  function  of  the  scalar  parameter  X,  which  weights 
the  previous  solution  with  *:he  current  optimal  solution  to  the  linear¬ 
ized  prob.em.  Denote  the  optimal  solution  to  (ii)  by  x*< 

-  „  £  , ,  ,  *.  £-1  ,  ,  *.£ 

Set  x  *  (1  -  X  )x  +  X  x  . 

4.  (Test  for  convergence.)  If  X*  =  0,  terminate  with  x5’1  as 
the  optimal  solution.  Otherwise,  set  x^  -  (1  -  X*)x^'^  +  X*x^  and  go 
to  step  5. 


-22- 


-l 

5.  Determine  tt  .  This  is  accomplished  as  in  step  3  of  •'he  cutting- 

plane  algorithm  by  sampling  from  the  distribution  of  5  and  solving  the 

resulting  linear  programs  parametrically  in  the  general  case,  or  expli- 
£ 

citly  calculating  tt  for  the  case  of  simple  recourse.  Go  to  step  2. 

Some  additional  remarks  about  the  algorithm  are  in  order.  The 
auxiliary  minimization  problem  that  is  required  in  step  3  requires  a 
considerable  amount  of  computation. 

An  alternative  convergence  criterion  is  available  at  step  A.  From 
results  in  [8]  an  upper  bound  on  the  difference  between  the  minimum 

4 

value  of  the  objective  function  z(x°)  and  the  current  value  z(x'“')  is 
given  by 

(12)  S(x£)  -  Ln  , 

where 

L  =  min  z(x^)  +  [c  -  n(x*1)T](x^  -  x^)  . 
h=l , . . . 

Thus  the  problem  may  be  terminated  when  (12)  is  less  than  some  pre¬ 
specified  constant. 

In  addition,  we  note  that  a  majority  of  the  constraints  in  the 
second-stage  problem  of  the  monthly  model,  i.e.  ,  those  of  the  form 

Za.  ,,  x.  £  a.  .x. .  , 
i-jk  Ljk  i  j  ij 

k^j 

and  some  of  the  first-stage  constraints  in  the  daily  model  are  of  a 
type  known  as  generalized  upper  bounds  [3].  These  are  constraints 
which  appear  schematically  as 


-23- 


That  is,  they  are  disjoint  sums  where  the  right-hand  side  is  a 

* 

nonnegative  constant.  A  special  algorithm  is  available  for  solving 
problems  of  this  form,  which  entails  carryi.-g  a  basis  where  size  is 
reduced  by  L,  the  number  of  generalized  upper-bound  constraints. 

Finally,  as  mentioned  earlier,  some  rounding  scheme  is  necessary 
to  convert  the  numbers  of  assigned  flights  to  integers.  One  scheme 
might  be  to  first  round  up  solutions  whose  fractional  parts  are  0.5 
or  greater  and  round  down  those  whose  fractional  parts  are  less  than 
0.5,  Then,  if  some  flying  hour  constraint  is  exceeded,  subtract  one 
from  a  sufficient  number  of  variables  to  satisfy  that  constraint.  The 
variables  to  be  reduced  should  be  chosen  from  among  those  previously 
rounded  up  and  may  be  those  whose  fractional  part  is  closest  to  0.5. 


* 

The  constraints  may  be  reduced  to  this  form  without  loss  of  gener¬ 
ality  by  dividing  each  coefficient  by  the  right-hand  side  and  rescaling 
the  variables. 
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V.  NUMERICAL  EXAMPLE 


We  conclude  with  a  small-scale  numerical  example  of  the  solution 
of  the  monthly  model.  The  problem  consists  of  two  aircraft  types  and 
two  routes.  The  data  were  chosen  for  convenience.  Solutions  based  on 
actual  data  will  be  reported  in  the  future. 

Ihe  data  for  the  problem  were  as  follows: 


Ai rcra  f t 

Flying  Hours 
(round  trip) 

Carrying  Capacity 
(in  tons) 

Cost  per  Flight 
(dollars) 

Type 

Route  1 

Route  2 

Route  1 

Route  2 

Route  1 

Route  2 

1 

24 

14 

50 

75 

1 

6000 

2 

49 

29 

20 

20 

4000 

Penalty  Cost 


($  per  ton) 

Route 

Excess  Demand 

Excess  Supply 

1 

500 

0 

2 

250 

0 

In  addition,  the  time  to  switch  between  either  route  was  5  hours 
for  aircraft  type  1  and  7  hours  for  aircraft  type  2.  The  incremental 
costs  for  these  moves  were  $1000  and  $1500  for  these  aircraft  types, 
respectively.  It  was  assumed  that  720  hours  of  flying  time  were  avail¬ 
able  for  each  aircraft  type.  This  can  correspond  to  3  planes  of  each 
type  with  -*n  average  utilization  of  8  hours  per  day. 

Finally,  demand  was  assumed  to  be  independently  lognormally  dis- 

•k 

tributed  on  each  route  with  the  following  parameters: 


The  current  programs  allow  only  for  the  lognormal  distribution. 
It  is  a  simple  matter  to  add  subroutines  to  generate  random  samples 
from  arbitrary  distributions  whose  cumulative  distribution  function 
can  be  tabulated.  The  lognormal  distribution  was  chosen  for  experi¬ 
mental  purposes  because  it  produces  nonnegattve  random  variables  which 
are  skewed  to  the  right  —  characteristics  actual  demand  patterns  are 
likely  to  exhibit. 
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Route 

U 

a 

1 

1000 

50 

2 

1500 

300 

The  following  solution  was  found  using  an  experimental  program  for 
the  cutting-plane  version  of  the  algorithm. 


Route 

Aircraft 

Type 

1 

2 

1 

16.5 

23.2 

2 

6.7 

0.0 

The  total  cost  of  the  entire  program  converged  after  4  terations 
(based  upon  a  termination  when  the  relative  decrease  in  costs  fell  below 
1  percent).  Each  second-stage  program  was  based  upon  drawing  a  sample 
of  25  observations  so  that  in  effect  about  100  linear  programming  prob¬ 
lems  were  solved  during  the  course  of  the  iterations.  Execution  time 
was  ten  seconds. 

For  comparison  purposes,  a  deterministic  linear  program  was  run 
with  demand  on  routes  assumed  to  be  equal  to  the  mean  used  in  the  log¬ 
normal  case.  This  yielded  the  following  results: 


Route 

Aircraft 

_  Type 

1 

2 

1 

18.3 

20.0 

2 

4.2 

0.0 

The  program  was  written  in  FORTRAN  and  employs  a  small,  all-in¬ 
core  linear  programming  subroutine.  Currently,  the  program  is  being 
modified  to  handle  much  larger  problems  and  will  be  described  in 
future  work.  A  version  of  the  algorithm  for  the  monthly  model  using 
the  gradient  algorithm  was  also  developed.  Preliminary  runs  indicate 
that  the  gradient  method  takes  considerably  longer  to  run  than  the 

cutting-plane  version  because  of  the  number  of  additional  programs  I 

that  must  be  solved  during  the  auxiliary  minimization  problem  described  ■ 

in  Sec.  IV.  I 


T 


It  should  be  noted  that  this  solution  differs  from  the  ca^e  where 
the  demands  are  stochastic.  It  might  be  supposed  that  the  stochastic 
'"-rsior.  would  always  assign  more  aircraft  as  a  hedge  against  uncer¬ 
tainty;  however,  this  is  not  the  case  even  in  this  simple  example  -- 
more  flig  ts  of  type  l  aircraft  to  route  1  are  made  in  the  determin¬ 
istic  version  than  in  the  stochastic. 
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